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Abstract 

\ 

The  measurement  of  liquid  water  in  snowpack  based 
on  the  electrical  path  length  difference  in  snowpack  depth 
in  the  region  of  high  dielectric  dispersion  of  water  is 
described.  The  measurement  system  consists  of  a  previously 
described  FM/CW  microwave  system  (Ellerbruck  and  Boyne, 

1980)  operating  at  two  frequency  bands  2-5  GHz  and  5-8 
GHz.  The  system  is  tested,  evaluated  and  itercompared 
with  a  dilution  measurement  technique  (Davis,  et.  al.  1985) 
in  both  a  laboratory  and  field  environment. 

It  is  found  that  the  two  methods  of  intercomparison 
are  equivalent  provided  the  snowpack  is  homogeneous  and 
free  of  stratigraphy.  The  accuracy  of  the  intercomparison 
is  ±2.5%  at  the  957«  level  of  confidence. 

The  intercomparison  of  the  methods  is  not  equivalent 
when  compared  in  an  inhomegeneous  stratified  snowpack. 

The  microwave  system  consistently  registers  a  higher  volumetric 
liquid  water  content  than  the  dilution  method.  This  discrepancy 
is  believed  due  to  differences  in  the  snowpack  volumes 
sampled.  In  particular,  it  is  believed  that  ponded  water 
above  impermeable  layers  regions  in  the  snowpack  is  responsible 
for  the  discrepancy.  These  regions  could  not  be  cared 
successfully  for  measurement  by  the  dilution  method. 
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MICROWAVE  DETERMINATION  OF  SNOWPACK  LIQUID  WATER  CONTENT 

Abstract 

The  measurement  of  liquid  water  in  snowpack  based  on  the  electrical 
path  length  difference  in  snowpack  depth  in  the  region  of  high  dielectric 
dispersion  of  water  is  described.  The  measurement  system  consists  of  a 
previously  described  FM/CW  microwave  system  (Ellerbruch  and  Boyne,  1980) 
operating  at  two  frequency  bands  2-5  GH z  and  5-8  GHz.  The  system  is 
tested,  evaluated  and  intercompared  with  a  dilution  measurement  technique 
(Davis,  et  al.  1985)  in  both  a  laboratory  and  field  environment. 

It  is  found  that  the  two  methods  of  intercomparison  are  equivalent 
provided  the  snowpack  is  homogeneous  and  free  of  stratigraphy.  The 
accuracy  of  the  intercomparison  is  ±  2.5%  at  the  95%  level  of  confidence. 

The  intercomparison  of  the  methods  is  not  equivalent  when  compared  in 
an  inhomogeneous  stratified  snowpack.  The  microwave  system  consistently 
registers  a  higher  volumetric  liquid  water  content  than  the  dilution 
method.  This  discrepancy  is  believed  due  to  differences  in  the  snowpack 
volumes  sampled.  In  particular,  it  is  believed  that  ponded  water  above 
impermeable  layered  regions  in  the  snowpack  is  responsible  for  the  discrep¬ 
ancy.  These  regions  could  not  be  cored  successfully  for  measurement  by  the 


dilution  method 


MICROWAVE  DETERMINATION  OF  SNCWPACK 


LIQUID  WATER  CONTENT 


I.  Introduction 

Research  on  the  interaction  of  electromagnetic  (ern)  waves  with  snow- 
pack  has  been  conducted  over  the  last  several  years  in  an  attempt  to  under¬ 
stand  the  response  of  em  waves  to  snowpack  and  to  use  this  interaction  to 
infer  or  predict  the  physical  state  of  snowpack,  such  as  depth,  water 
equivalence,  and  liquid  water  content. 

Snowpack  is  a  dynamic  medium  continually  undergoing  changes  in  dens¬ 
ity,  grain  size  and  shape,  stratigraphy,  water  equivalence,  and  liquid 
water  content.  These  changes  are  driven  by  radiation  and  advection  proces¬ 
ses  and  storm  patterns.  Within  this  spectrum  of  changing  snow  conditions, 
the  physical  characteristics  of  liquid  water  content  and  density  have  con¬ 
siderable  impact  on  the  interaction  of  electromagnetic  waves  with  snow¬ 
pack.  The  modification  of  the  electromagnetic  interaction  is  brought  about 
by  the  differences  in  the  dielectric  properties  of  water  and  ice.  At  a 
frequency  of  10  GHz,  the  relative  dielectric  constants  of  water  and  ice  at 
0°C  are  40  and  3.15,  respectively,  and  the  rates  of  attenuation  are  4  x  103 
Np/meter  and  5  x  10~2  Np/meter,  respectively  (Stiles  and  Ulaby,  1980a). 

The  relative  dielectric  constant  of  dry  snow  having  a  density  of  350  kgm/m 
is  1.8.  Depending  on  the  liquid  water  content,  snowpack  can  change  from  a 
relatively  lossless  dielectric  to  a  high  loss  dielectric.  The  changes  in 
the  dielectric  properties  can  affect  the  response  of  both  active  and 
passive  microwave  systems. 
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Several  experiments  have  been  conducted  to  study  the  interaction  of 
electromagnetic  waves  with  snowpack  (Ellerbruch  and  Boyne,  1980;  Stiles  and 
Ulaby,  1980a, b;  Ulaby  and  Stiles,  1980;  Chang  et  al. ,  1980;  Matzler  et  al. , 
1980;  Linlor,  et  al.  1980).  Active  radar  measurements  have  been  reported 
at  frequencies  in  the  range  1  to  37  GHz,  while  passive  radiometric  measure¬ 
ments  have  been  made  in  the  range  5  to  94  GHz.  These  experiments  have 
shown  that  the  interaction  of  em  waves  with  snowpack  is  governed  by  both 
surface  and  volumetric  effects  of  the  snowpack  and  may  be  influenced  by  the 
underlying  ground  surface.  The  degree  to  which  volumetric  effects  play  a 
role  is  dependent  on  the  snow  wetness,  stratigraphic  layering,  grain  size, 
and  the  observing  frequency. 

The  major  difficulty  with  these  experiments  is  that  the  snowpack  scene 
is  too  con^lex  to  allow  a  rigorous  quantitative  evaluation  of  the  electro¬ 
magnetic  response  to  snowpack  water  equivalence  and  liquid  water  content. 

In  all  experiments  the  snow  medium  was  Inhomogeneous  and  stratified.  It 
contained  snow  in  several  different  stages  of  metamorphism  ranging  from 
depth  hoar  near  the  ground-snow  interface  to  melt-freeze  near  the  snow-air 
interface. 

The  generation  of  liquid  water  in  snowpack  is  dynamic  and  commonly 
exhibits  spatial,  volumetric  and  temporal  variations,  especially  in  layered 
snowpack  (Colbeck,  1978;  Wankiewicz,  1979).  In  situ  measurements  are  dif¬ 
ficult  and  the  permeability,  thickness  and  areal  extent  of  buried  crusts 
and  ice  layers  is  unknown.  Thus,  while  experiments  in  wet  snowpack  have 
been  conducted,  the  inability  to  establish  spatial  and  temporal  liquid 
water  profiles  in  the  pack  has  limited  the  interpretation  of  the  results 
(Kong,  et  al.  1980). 
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Experiments  specifically  designed  to  measure  the  dielectric  properties 
of  snow  containing  liquid  water  under  controlled  laboratory  conditions  have 
been  conducted  (Hallikainen  et  al.  ,  1985;  Sweeny  and  Colbeck,  1974). 
However,  the  measurements  were  made  to  compare  the  dielectric  properties  of 
wet  snow  with  various  dielectric  mixing  models  rather  than  developing 
methodology  to  measure  liquid  water  content  in  a  snowpack  environment. 
Direct  measurements  of  liquid  water  content  in  snowpack,  which  are  con¬ 
sidered  reliable,  are  the  freezing  calorimeter  technique  (Jones  et  al. , 
1980),  the  capacitance  technique  (Ambach  and  Denoth,  1972,  1980),  and  the 
dilution  method  (Davis  et  al. ,  1985).  These  techniques  sample  only  a 
limited  volume  of  snowpack,  which  is  not  necessarily  representative  of  the 
entire  snowpack. 

The  measurements  described  in  this  report  represent  attest  to  deter¬ 
mine  liquid  water  content  of  a  large  snowpack  volume  and  account  for  varia¬ 
tions  of  liquid  water  within  the  snowpack.  A  measurement  system  is 
described  which  is  sensitive  to  the  electrical  path  length  between  the  top 
and  bottom  of  the  snowpack.  A  three  conq>onent  model  for  snowpack  is 
developed  where  the  electrical  path  length  through  the  snowpack  is  differ¬ 
ent  in  two  different  frequency  bands  due  to  the  dielectric  dispersive 
property  of  water.  These  differences  are  detectable  and  form  the  basis  of 
a  measurement  instrument  to  determine  liquid  water  content  in  snowpack. 

The  purpose  of  the  study  is  to  test  and  evaluate  the  model  by  develop¬ 
ing  a  microwave  system  to  measure  liquid  water  in  snowpack.  The  microwave 
measurements  are  intercompared  with  a  dilution  measurement  method  which  has 
been  described,  tested  and  analyzed  previously  (Davis,  et  al.  1985). 
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II.  Electrical  Path  Length  Model  (EPLM) 

Consider  a  snowpack  consisting  of  ice,  air,  and  liquid  water.  The 
volumetric  liquid  water  content  is  restricted  to  be  not  greater  than  10%. 
This  regime  is  typical  of  a  well-drained  ‘'ripe”  snowpack. 

The  dielectric  properties  of  the  snowpack  are  treated  as  a  lumped 
homogeneous  system.  That  is,  the  three  components  of  the  mixture  (ice, 
air,  and  water)  are  lumped  into  layers  each  of  which  is  considered  dielec¬ 
trically  homogeneous  (Colbeck,  1980).  The  porosity  is  defined  as  <J>  *  9a 
+  0 ,  where  0a  and  8  are  the  volumentric  contents  of  air  and  water, 
respectively,  and 

p  *  (1  -  $)  p  +  0p  +  9  p  »  (1  -  <p )  p  +  0p 
s  i  w  a  a  iw 

where 

pg  *  density  of  snow 

p.^  ■  density  of  ice  (0.917  Mg/m3) 

pw  *  density  of  water 

pa  *  density  of  air 

The  relative  dielectric  constants  of  snow,  ice,  air,  and  water  are 
eg,  e^,  and  ew,  respectively,  and  the  electrical  path  length  through 
a  vertial  section  of  the  snowpack  is 


/  e  d 
s  s 


/e  d  +  /e .d.  +  /e  d 
a  a  1  i  w  w 


(1) 


For  a  unit  cross  sectional  area,  the  depths  of  each  constituent  are 

0. 

t 

dy  *  9  ■  depth  of  water;  and  dj_  ■  1  -  $  is  the  depth  of  ice.  Since  1  - 


numerically  equal  to  their  volumetric  contents;  da  *  9a  ■  depth  of  air; 


+  9  +  9a  ■  1,  dg  represents  unit  snow  depth. 


Thus, 
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(2) 


/e  *  /e  0  +  /e  (1  -<(>)+  /e  9 

s  a  a  i  w 

In  the  microwave  region  of  1-12  GHz,  ea  and  are  frequency  independ- 
I  ent.  ew  is,  however,  frequency  dependent. 

Equation  (2)  is  the  basis  of  comparison  with  other  data.  Figure  1  is 
a  comparison  with  the  data  of  Cumming  (1952)  and  Sweeny  and  Colbeck  (1974) 
for  dry  snow  (see  Colbeck,  1980).  Figure  2  compares  the  data  of  Ambach  and 
Denoth,  and  the  electrical  path  length  model  for  pg  =*  0.32  Mg/m3,  ew  = 

87.91  corresponding  to  a  frequency  of  20  MHz,  which  was  used  by  Ambach  and 
Denoth.  Figure  3  is  a  comparison  of  the  Polder-Van  Santen  model  calcula¬ 
tion  (0-10%)  developed  by  Colbeck  (1980)  and  the  electrical  path  length 
model  for  pg  *  0.32  Mg/cm3  ,  with  ew  *  87.91.  The  agreement  between  the 
electrical  path  length  model  and  the  Ambach  and  Denoth  data  is  quite 
reasonable.  Similarly,  the  comparison  of  Figure  3  shows  a  maximum  discrep¬ 
ancy  of  6%  at  0  *  10%. 

III.  Experimental  Procedure 

The  FM-CW  radar  operating  in  the  range  1-12  GHz  is  sensitive  to  the 
electrical  path  length  in  a  snowpack.  This  observation  has  been  used  to 
measure  snowpack  water  equivalence  of  dry  snowpack  by  measuring  the  fre¬ 
quency  separation  between  reflections  of  swept  frequency  microwaves  at  the 
top  and  bottom  of  the  snowpack  (Ellerbruch,  and  Boyne  1980).  The  governing 
equation  is  restated  here 

Af  *  K/F  d  (3) 

s  s 

where  ea  is  the  dielectric  constant  of  sncw,  dQ  is  the  depth  of  snow 
4(f2-fx)fn 

and  K  -  - — -  is  a  system  constant  with  f2“fi  the  microwave  frequency 
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band  actually  swept,  fn  the  sweep  frequency  and  c  the  speed  of  light. 
Writing  equation  (1)  for  two  frequency  bands,  a  and  8,  we  have 

/e  d  */e  d  +  /e.d.  +  /e  d 
sa  s  a  a  i  i  ua  u 


/e  .  d  =  /e  d  +  /e.d,  +  /e  _d 
s8  s  a  a  i  i  uP  in 


Subtracting  these  equations,  we  have 


f/e  -  /e  a)d 
^  _  2  sa _ s8;  s 

U)  /  “ 

/e  -  ve  Q 
aj  a  0)8 

ea)a  and  e^g  are  determined  by  averaging  the  dielectric  constant  over 
the  frequency  band  of  interest,  i.e. , 

h 

o)“  c  df 

-  i.  *  +  «^„>2 


(4) 


e  -  e 

£  (f7  ~  f.)  -  £  +  - 7— 

W  1  1  0)»  f2  f 


(5) 


where  -87.91,  *  4.9,  and  fm  *  8.51  x  lCr  Hz  at  a  temperature 

of  0°C. 

Substituting  equation  (3)  into  equation  (4)  gives 

.  4VK«  -  -VS 

a  *  ~  - —  - 

/ £  -  /£  „ 

awt  0)8 

where  Afa  and  Afg  are  measured  on  the  same  volume  of  snowpack. 

Assuming  the  same  system  constant  K  for  both  frequency  bands,  a 


(6) 


frequency  difference  of  30  to  60  Hz  can  be  expected  for  a  1%  volume  of 
liquid  water  assuming  the  two  microwave  frequency  bands  to  be  2-5  GHz,  and 
5-8  GHz,  the  sweep  frequency  150  Hz  and  a  snow  depth  0.5  m  to  1  m. 


IV  Estimation  of  Measurement  Accurac 


Equation  6  can  be  analyzed  to  estimate  the  accuracy  of  the  microwave 
measurement  method.  The  error  associated  with  the  determination  of  d^  is 

H,  ,  6K  .  2  li(Afae) 

d  K  fAf  ~  Af 

u>  !  a  6 

where  it  has  been  assumed  that  K^  =  Kg  and  that  the  errors  associated 
with  Afa  and  Afg  are  comparable.  The  error  in  K  is 

I  5K  6  ( f  7  ~f ,  )  . 

K  "  f2-f,  f 

Measurements  of  the  voltage  controlled  frequency  sweep,  f2  “  f  j. .  are  repro 
ducible  in  the  laboratory  to  an  accuracy  of  ±  0.1%  and  fn  is  reproducible 
to  ±  0.1%  (99%  confidence  interval).  Therefore, 

J 6K/Kj  -  0.2%  . 

Laboratory  reproducibility  of  Afa  or  Afg  is  ±  5  Hz.  For  Afa  * 

30  Hz  per  1%  change  in  the  liquid  water  content,  we  see  that  the  error 
associated  with  6K/K  is  negligible  in  comparison  to  | 6 (Af ) //Af a-Af g | 


6  d 

6Af 

0) 

d 

0) 

*  L 

Af  - 

0.  P 

Figure  4  shows  the  relative  error  as  a  function  of  liquid  water  content, 
the  error  being  ±  33%  at  d^/dg  =  1%  and  3.3%  at  d^/dg  =  10%  for  the 
micrcwave  frequency  bands  of  2-5  GHz  and  5-8  GHz. 


Dilution  Method: 


The  analytical  estimate  of  errors  with  the  dilution  method  has  been 
analyzed  by  Davis  et  al.  (Davis,  et  al.  1985),  They  estimate  the  errors 
due  to  the  combined  measurements  of  the  conductivity  and  mass  of  the  solu¬ 
tion,  snow  and  the  mixture  to  be 

X2 

|«  |  ■  10'2  S/M  +  1.3  x  10-2  Xu  +  6.5  x  10-  J7J 

where  S  is  the  mass  of  stock  solution,  M  is  the  mass  of  the  snow  sample, 
and  X^  is  the  liquid  water  mass  fraction.  For  an  S/M  ratio  of  between 
0.4  and  1.0,  the  results  of  72  laboratory  tests  give  a  mean  absolute  error 
of  1.2%  for  Xy. 

V.  Description  of  Experiments 

Experiments  were  conducted  in  two  phases:  laboratory  cold  room 
measurements  and  field  measurements.  The  measurements  were  designed  to 
verify  the  measurements  of  liquid-water  content  in  snowpack  with  the  two 
frequency  band  FM-CW  system. 

A.  Laboratory  Experiments 

Snow  samples  were  prepared  in  an  insulated  box  of  dimensions  1  x  1  x 
0.5  m  (Fig.  5).  The  box  had  a  false  bottom  which  was  composed  of  ceramic 
porous  plates  having  an  apparent  porosity  of  54%.  The  ceramic  plates  were 
covered  with  a  metal  screen.  The  bottom  of  the  box  contained  an  outlet  to 
drain  off  excess  water  in  the  snowpack  and  to  provide  tension  on  the  snow- 
pack  to  drain  the  pack  to  a  specified  volume  of  liquid  water.  The  box  was 
placed  in  a  laboratory  cold  room  capable  of  temperature  regulation  from  +10 


to  -18°C 


Microwave  antennas  were  situated  to  reflect  the  radiation  into  the 


box.  The  radiation  intensity  in  the  box  was  measured  to  determine  its  dis¬ 
tribution  (Fig.  6).  Natural  snow  was  collected  and  sifted  into  the  box 
through  a  screen  whose  average  aperture  was  5  mm.  The  snow  was  then 
covered  with  styrofoam  and  kept  at  constant  temperature  of  -1  to  -56C  for 
several  days  to  allow  the  snow  to  settle  and  stabilize.  Liquid  water  was 
generated  by  bringing  the  temperature  to  0°C  and  exposing  an  aluminum  plate 
placed  on  top  of  the  pack  to  infrared  radiation  generated  from  one  or  two 
lamps.  The  length  of  time  the  snow  was  irradiated  was  determined  by  the 
volumetric  water  content  desired.  Several  tests  of  this  process  were  made 
to  determine  the  degree  of  uniformity  of  liquid  water  content  within  the 
box  (Fig.  7). 

Microwave  measurements  were  made  on  the  snow  sample.  The  response  of 
the  FM/CW  signal  was  displayed  on  a  spectrum  analyzer  and  subsequently 
recorded  on  an  x-y  plotter.  Figure  8  shows  the  response  of  both  the  2-5 
GHz  and  5-8  GHz  signals  for  a  typical  run. 

Immediately  afterwards,  three  0.054  m  diameter  core  samples  ranging 
from  0.4  m  to  0.5  m  in  depth  were  taken  for  direct  measurement  of  liquid 
water  by  the  dilution  method.  The  core  samples  were  distributed  in  such  a 
way  that  three  sets  of  measurements  could  be  made  before  the  snow  had  to  be 
redistributed  in  the  box. 

B.  Field  Experiments 

Field  experiments  were  conducted  at  the  summit  of  Cameron  Pass  at  an 
elevation  of  3200  m.  Measurements  were  conducted  in  the  late  Spring  of 
1984  and  1985.  The  sncwpack  was  at  or  near  equitemperature  and  the  experi¬ 
ments  were  conducted  in  a  large  meadow  on  level  terrain.  The  general  pro- 


cedure  for  the  field  experiments  was  to  dig  an  observation  pit,  generally 
1.5  to  2  m  deep,  so  that  density,  stratigraphy,  grain  size  and  type  and 
liquid  water  content  could  be  measured  (Fig.  9).  The  microwave  system  was 
set  up  near  the  south  wall  of  the  pit.  In  all  field  measurements  the 
bottom  of  the  snowpack  sample  was  determined  by  the  insertion  of  a  metal 
plate  from  which  the  microwave  signal  was  reflected.. 

Microwave  measurements  were  made  on  a  1  m2  cross  section  of  snowpack. 
Since  the  top  0.5  m  of  snowpack  had  considerable  structure,  including  ice 
lenses  (Fig.  9),  measurements  were  also  made  on  the  snowpack  after  removing 
the  top  0.5  m  to  0.7  m  of  snow.  The  pack  at  this  depth  was  generally 
uniform  and  free  of  stratigraphy  for  the  remainder  of  the  day  of  observa¬ 
tion.  However,  night  radiation  loss  usually  froze  the  exposed  surface  such 
that  the  generated  liquid  water  was  spatially  nonuniform  on  the  day  follow¬ 
ing  exposure. 

Direct  measurements  of  the  liquid  water  content  were  made  using  the 
dilution  technique.  Three  to  five  core  samples  0.054  m  in  diameter  and  0.4 
m  to  0.5  m  in  depth  were  taken  immediately  after  each  micrtwave 
measurement.  After  each  microwave  and  dilution  measurement  comparison  was 
completed,  a  new  sample  volume  of  snow  was  prepared. 

A  pit  profile  recorded  stratigraphy,  grain  size  and  type,  density  and 
temperature  at  various  depths  (Fig.  10  and  11).  Density  measurements  were 
made  with  a  glacial  sampler  after  each  measurement  rather  than  relying  on 
the  core  measurement  since  the  diameter  of  the  core  sampler  was  small 
enough  to  induce  systematic  errors  in  the  density  (Fames  et  al. ,  1982). 


The  data  for  both  the  laboratory  and  field  experiments  which  compare 
the  dilution  and  microwave  measurements  of  liquid  water  content  in  snowpack 
can  be  found  in  the  Appendix.  Tables  1,  2  and  3  summarize  the  laboratory 
and  field  data.  Column  3  of  the  tables  represents  the  liquid  water  content 
determined  by  the  dilution  method  while  column  4  is  the  liquid  water  con¬ 
tent  determined  by  the  microwave  method  on  the  same  snowpack  sample. 

Column  5  represents  the  deviation  and  column  6  the  variance  of  each  inter- 
comparison.  The  data  in  Tables  1  and  2  show  reasonable  agreement  between 
the  two  measurement  methods.  In  these  comparisons,  emphasis  was  placed  on 
sampling  well-drained  uniform  snow  which  was  free  of  stratigraphic  layering 
and  had  a  reasonably  uniform  grain  structure  typical  of  a  moist  equitem- 
perature  snowpack. 

The  comparisons  of  Table  3,  however,  were  taken  in  snow  that  had 
considerable  stratigraphy  and  large  variations  in  grain  structure  between 
layers.  It  is  evident  that  the  two  measurement  methods  are  not  comparable 
in  this  situation.  Figures  10  and  11  show  the  snowpit  data  for  runs  1 
through  6  of  Table  3.  The  numbers  to  the  right  of  the  figure  correspond  to 
the  first  six  runs  made.  The  depth  of  snow  sampled  for  each  run  is  also 
shown.  Air  temperature  was  4°C  to  5°C  and  surface  melting  was  in  prog¬ 
ress.  However,  the  stratigraphy  in  the  pack  maintained  its  integrity 
throughout  the  day.  Nighttime  temperatures  were  below  freezing  creating  a 
hard  frozen  crust  on  the  sncw  surface  and  occasional  surface  hoar.  Because 
of  the  ice  lenses  and  sun  crust  it  was  difficult  to  obtain  core  samples  for 
dilution  measurements.  Sampling  could  be  done  only  where  the  stratigraphy 
was  weakened  by  percolation  of  water.  Therefore,  while  the  microwave 
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VI.  Experiment  Results 


The  data  for  both  the  laboratory  and  field  experiments  which  compare 
the  dilation  and  microwave  measurements  of  liquid  water  content  in  snowpack  • 
can  be  found  in  the  Appendix.  Tables  1,  2  and  3  summarize  the  laboratory 
and  field  data.  Column  3  of  the  tables  represents  the  liquid  water  content 
determined  by  the  dilution  method  while  column  4  is  the  liquid  water  con¬ 
tent  determined  by  the  microwave  method  on  the  same  snowpack  sample. 

Column  5  represents  the  deviation  and  column  6  the  variance  of  each  inter¬ 
comparison.  The  data  in  Tables  1  and  2  show  reasonable  agreement  between 
the  two  measurement  methods.  In  these  comparisons,  emphasis  was  placed  on 
sampling  well-drained  uniform  snow  which  was  free  of  stratigraphic  layering 
and  had  a  reasonably  uniform  grain  structure  typical  of  a  moist  equitem- 
perature  snowpack. 

The  comparisons  of  Table  3,  however,  were  taken  in  snow  that  had 
considerable  oCratigraphy  and  large  variations  in  grain  structure  between 
layers.  It  Is  evident  that  the  two  measurement  methods  are  not  comparable 
in  this  situation.  Figures  10  and  11  show  the  snowpit  data  for  runs  1 
through  6  of  Table  3.  The  numbers  to  the  right  of  the  figure  correspond  to 
the  first  six  runs  made.  The  depth  of  snow  sampled  for  each  run  is  also 
shown.  Air  temperature  was  4°C  to  5°C  and  surface  melting  was  in  prog¬ 
ress.  However,  the  stratigraphy  in  the  pack  maintained  its  integrity 
throughout  the  day.  Nighttime  temperatures  were  below  freezing  creating  a 
hard  frozen  crust  on  the  snow  surface  and  occasional  surface  hoar.  Because 
of  the  ice  lenses  and  sun  crust  it  was  difficult  to  obtain  core  samples  for 
dilution  measurements.  Sampling  could  be  done  only  where  the  stratigraphy 
was  weakened  by  percolation  of  water.  Therefore,  while  the  microwave 
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measurements  sampled  the  entire  volume  Including  water  which  was  ponded 
above  the  stratigraphic  layers,  the  dilution  samples  were  taken  in  a  part 
of  the  volume  where  the  pack  was  draining. 

The  question  of  equivalence  of  these  two  methods  of  measurement 
involves  an  analysis  of  variance  of  the  paired  comparisons.  Tables  4,  5 
and  6  summarize  the  tests  of  significance  using  the  F  test  statistic 
(Mandel,  J. ,  1964).  The  statistical  test  compares  two  models,  one  in  which 
Xy(FM)  +  e  *  XV(DL)+  5  where  e  and  6  are  the  random  errors  associated 
with  each  measurement  method,  and  one  in  which  XyCFM)  +  e  *  K  +  X^DL) 

+  6  where  K  represents  a  constant  difference  between  the  results.  The 
hypothesis  to  be  tested  is  that  K  ■  0. 

Table  4,  5  and  6  present  the  analyses  for  the  data  presented  in  Tables 
1,  2  and  3.  The  analyses  show  that  the  test  for  K  *  0  is  not  significant 
at  the  5%  level  for  the  data  presented  in  Tables  4  and  5,  but  is  signifi¬ 
cant  for  the  data  presented  in  Table  6.  Therefore,  K  -  0  for  the  data  in 
Tables  4  and  5  and  the  measurement  of  liquid  water  by  the  microwave  and 
dilution  methods  measurements  are  equivalent.  However,  K  *  0  for  the  data 
in  table  6  and  the  paired  comparison  sh<xr  a  definite  bias.  The  average 
variance  for  the  combined  paired  comparisons  of  Table  4  and  5  is  1.61  and 
the  standard  deviation  1.27.  Thus  the  uncertainty  associated  with  the 
comparison  at  the  95Z  confidence  level  is  ±  1.96o  or  ±  2.5%. 

VII  Discussion; 

It  is  apparent  that  the  ability  to  compare  liquid  water  measurement 
methods  will  depend  to  a  large  extent  on  the  character  of  the  snowpack.  If 
the  pack  is  homogeneous  and  free  of  stratigraphy,  reasonable  success  of 
intercomparisons  can  be  achieved.  If,  however,  the  snowpack  is  inhomogene- 


ous  with  several  stratigraphic  layers,  it  may  be  difficult  to  achieve  even 
order  of  magnitude  accuracy  due  to  variations  snow  texture,  liquid  water 
content,  liquid  water  distribution  and  differences  in  the  sample  volumes 
measured  by  the  different  methods.  Difficulties  of  this  type  have  been 
reported  in  intercomparisons  of  several  electrical  sensors  operating  at 
frequencies  ranging  from  1  MHz  to  1.3  GHz  (Denoth,  et  al.  1984).  These 
senors  showed  agreement  when  intercompared  in  wet  sand  and  prepared  homo¬ 
geneous  wet  snow  to  within  ±  1%,  but  discrepancies  of  up  to  10%  were 
reported  in  natural  inhomogeneous  snowpack. 

Water  flow  through  layered  snowpack  media  is  understood  only  qualita¬ 
tively  and  in  situ  measurements  of  liquid  water  distribution  within  the 
snow  are  difficult.  Wankiewiez  (1979)  has  proposed  that  flow  through 
layered  snow  be  modeled  in  such  a  way  that  the  boundary  between  two  snow 
layers  of  different  textures  and  densities  will  impede,  have  no  effect  on 
or  accelerate  downward  flow  depending  on  the  diffential  effect  of  pressure 
Immediately  above  and  below  the  boundary.  Using  this  model,  he  discusses 
the  effect  of  ice  layers  on  flow,  treating  them  as  confound  horizons. 
Various  likely  flow  patterns  resulting  from  different  combinations  of  fine 
and  coarse-grain  snow  are  presented,  all  of  which  result  in  spatially 
Inhomogeneous  distributions  of  meltwater  in  snowpack.  A  limited  number  of 
observations  (Wakahama,  1963,  1968,  Wankiewiez  and  de  Vries,  1978)  tend  to 
support  at  least  qualitatively,  the  proposed  model.  It  is,  therefore,  not 
surprising  that  the  two  methods  of  measuring  liquid  water  in  stratified 
snowpack  do  not  agree.  One  can  expect  ponding  of  liquid  water  above 
impermeable  layers  which  could  not  be  examined  with  the  core  sampling 
technique  used  for  dilution  measurements. 


Several  intercoqparisons  between  the  microwave  and  freezing  calorimet¬ 
er  methods  of  measuring  liquid  water  content  were  attempted.  The  freezing 
calorimeter  method  was  initially  intended  as  the  standard  for  the  inter¬ 
comparisons  since  the  basic  principles  of  the  method  are  well  understood. 
However,  the  technique  is  difficult  to  use  unless  one  is  intercomparing 
single  0.2  Kg  to  0.5  Kg  samples.  The  measurement  of  a  single  sample  with 
this  technique  is  lengthy  and  several  samples  (3  to  5)  are  needed  to  obtain 
a  liquid  water  profile.  Also,  the  technique  does  not  always  give  success¬ 
ful  results.  Special  precautions  must  be  taken  to  ensure  complete  mixing 
of  the  working  fluid  and  the  snow  sample.  Finally,  the  liquid  water  dis¬ 
tribution  in  a  reasonably  homogeneous  snowpack  may  vary  significantly, 
requiring  several  profiles  to  be  taken  to  account  for  these  variations. 

O 

Figure  12  shows  the  results  of  a  prepared  smoothed  1.5  m  surface  which  was 
sprayed  with  a  fine  mist  of  dye  and  left  exposed  to  solar  radiation  for 
approximately  10  minutes.  The  dark  splotches  are  where  meltwater  is  being 
generated.  It  is  unlikely  that  a  snow  sample  of  0.2  to  0.5  Kg  of  snow  used 
in  the  freezing  calorimeter  method  would  adequately  represent  even  the 
surface  meltwater  content. 

Initially  we  proposed  to  use  two  sets  of  frequency  bands  (2-5  and  5-8 
GHz,  and  2-6  and  8-12  GHz)  for  the  microwave  measurement.  However,  we  were 
not  successful  in  obtaining  adequate  response  from  the  8-12  GHz  system 
except  in  dry  snow.  Resent  results  on  the  dielectric  loss  in  this  frequen¬ 
cy  range  as  a  function  of  snow  wetness  (Hallikainen,  et  al. ;  1985)  indicate 
Chat  the  attenuation  at  10  GHz  is  approximately  42  db/m  per  1%  liquid  water 
content  (Hallikainen  et  al.  ,  Fig.  15). 


VIII  CONCLUSIONS 


The  objective  of  this  study  was  to  test  and  evaluate  the  performance 
of  a  swept  frequency  microwave  system  in  measuring  the  liquid  water  content 
of  snowpack,  and  to  intercompare  this  system  with  other  measurement  methods 
in  both  a  laboratory  and  field  environment. 

The  study  concludes  that  liquid  water  content  in  snowpack  can  be 
determined  using  the  microwave  system  in  the  2-5  GHz  and  5-8  GHz  frequency 
bands.  An  analysis  of  the  errors  indicates  that  the  measurement  can  be 
made  with  an  accuracy  of  ±10%  for  volumetric  liquid  water  content  greater 
than  2.5%  to  3.5%  depending  on  the  choice  of  system  parameters. 

Intercom) arison  of  the  microwave  and  dilution  method  of  measurement 
were  conducted.  It  was  found  that  the  two  methods  were  equivalent  provided 
the  snowpack  was  homogeneous  and  free  of  stratigraphy.  The  mean  standard 
deviation  of  the  intercomparison  for  the  laboratory  and  field  measurements, 
combined  was  1.27%,  thus,  the  error  associated  with  the  intercomparison  was 
±  2.5%  at  the  95%  level  of  confidence. 

It  was  found  that  the  intercomparison  was  not  equivalent  in  a  strati¬ 
fied  inhomogeneous  snowpack  and  that  a  definite  bias  was  evident,  the 
microwave  measurement  being  consistently  higher  than  the  dilution  measure¬ 
ment.  The  bias  is  attributed  to  the  difference  in  the  sample  volumes  of 
the  two  measurements,  the  microwave  measurement  sensing  a  1x1x0. 5m 
volume  and  the  dilution  measurement  sampling  three  cores  0.054  m  in  diamet¬ 
er  and  0.5  m  depth.  Furthermore  core  samples  could  not  be  acquired  where 
the  stratigraphic  layers  were  impermeable.  These  impermeable  regions  are 
where  ponding  of  liquid  water  is  most  likely  to  occur. 


The  attenuation  of  the  microwave  signal  in  the  region  8-12  GHz  was 
found  to  be  too  large  to  identify  the  bottom  of  the  snowpack  sample  when 
liquid  water  was  present.  The  attenuation  was  not  measured  with  great 
precision,  but  is  consistent  with  the  attenuation  reported  by  Hallikainen 
(Hallikainen  et  al.  1985). 

The  inter comparison  of  measurements  with  the  freezing  calorimeter  was 
abandoned  since  the  time  involved  in  obtaining  a  representative  liquid 
water  profile  in  the  snowpack  sample  volume  was  found  to  be  prohibitive. 
This  method  is  useful  if  small  volumes  are  intercompared  under  well 
controlled  conditions,  but  loses  its  utility  when  the  measurement  time, 
such  as  is  necessaiy  for  liquid  water  profile  measurements,  becomes 
comparable  to  the  time  in  which  snow  conditions  such  as  generation  of 
liquid  water  and  breakdown  of  layers  can  be  expected  to  change. 

Observation  of  surface  melt  using  dye  indicates  that  the  melt  process 
is  not  uniform  and  may  depend  on  the  past  micrometeorological  conditions 
such  as  eddy  currents  resulting  in  differences  in  the  surface  snow  texture 
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Table  1.  Liquid  water  in  snowpack  microwave-dilution  comparison 


Microwave  Frequency  Bands 
2-5  GHz  and  5-8  GHz 

Laboratory  Data 


Homogeneous  Snowpack. 


Date 

X  (DL) 
(vol  %) 

X  (FM) 

(vol  %) 

Deviation 

Variance 

1 

5  Feb 

2.96 

3.68 

-0.72 

0.518 

2 

6  Feb 

3.33 

5.56 

-2.23 

4.973 

3 

7  Feb 

1.26 

0.18 

1.08 

1.166 

4 

15  Mar 

3. 10 

3.22 

-0.12 

0.014 

5 

15  Mar 

1.25 

3.36 

-2.11 

4.452 

6 

20  Mar 

1.73 

2.10 

-0.37 

0.137 

7 

21  Mar 

4.10 

3.21 

0.89 

0.792 

8 

21  Mar 

5.95 

3.31 

2.64 

6.970 

9 

10  Apr 

1 .64 

1.21 

0.43 

0.185 

I  =  -0.51  19.207 


Mean  a  =  1.461 


Table  2.  Liquid  water  in  snowpack 
microwave-dilution  comparison. 

Microwave  Frequency  Bands 
2-5  GHz  and  5-8  GHz 

Field  Data 

X  (DL)  X  (FM) 

(vol  %)  (vol  X) 


Date 


Deviation 


Variance 


v000^40>Ln4>U>N) 


Table  3.  Liquid  water  in  sncwpack  microwave -dilution  comparison. 

Micrcwave  Frequency  Bands 
2-5  GHz  and  5-8  GHz 

Field  Data 


Date 

X  (DL) 
(vol  X) 

X  (FM) 
(vol  X) 

Deviation 

Variance 

24  Apr 

1.36 

10.0 

-8.64 

Top 

of  pack 

24  Apr 

1.24 

4.19 

-2.95 

Top 

15  cm  removed 

24  Apr 

1.81 

5.60 

-3.79 

Top 

cm  removed  25 

1  May 

3.73 

29.5 

-25.77 

Top 

of  pack 

1  May 

2.58 

13 

-10.42 

Top 

10  cm  removed 

1  May 

4.64 

20 

-15.36 

Top 

cm  removed  40 

7  May 

2.76 

8.70 

-5.94 

Top 

of  pack 

8  May 

3.62 

7.4 

-3.78 

Top 

of  pack 

8  May 

2.68 

10.3 

-7.62 

Top 

of  pack 

10  May 

3.37 

6.6 

-3.23 

Top 

cm  removed  30 

l  -  87.52 


21 


Table  4.  Liquid  water  in  snowpack  microwave-dilution  comparison, 

Analysis  of  Variance 
Laboratory  Data 


Homogeneous  Snowpack 


X  (DL) 
Date  (vol  X) 


1 

5 

Feb 

2.96 

2 

6 

Feb 

3.33 

3 

7 

Feb 

1.26 

4 

15 

Mar 

3.10 

5 

15 

Mar 

1.25 

6 

20 

Mar 

1.73 

7 

21 

Mar 

4.10 

8 

21 

Mar 

5.95 

9 

10 

Apr 

1.64 

X  (FM) 
(vol  %) 

Deviation 

V 

3.68 

-0.72 

5.56 

-2.23 

0.18 

1.08 

3.22 

-0.12 

3.36 

-2.11 

2.10 

-0.37 

3.21 

0.89 

3.31 

2.64 

1.21 

0.43 

y 

=  -0.51 

Deviation 


V2 

v-v 

(v-v)2 

0.5184 

-0.6633 

0.44001 

4.9729 

-2.173 

4.7234 

1.1664 

1.1367 

1.2920 

0.0144 

-0.0633 

0.00401 

4.4521 

-2.053 

4.2162 

0.1369 

-0.3133 

0.09818 

0.7921 

0.94667 

0.89618 

6.9696 

2.6967 

7.2720 

0.1849 

0.48667 

0.23684 

0.00 

Model 

Residual 

DF 

Sum  of  Squares 
of  residuals 

Mean 

square 

v  (FM) 

V 

-  X  (DL)  - 

V 

0  +  (e-5) 

9 

OO 

o 

CN 

• 

o> 

f—4 

n 

CM 

> 

1>-J 

2.134 

X  (FM) 

V 

-  X  (DL)  - 

V 

K  +  (e-6) 

8 

l  (v-v)2  -  19.179 

2.397 

Difference 

1 

£a  -  0.029 

0.029 

F  -  0.021 


Ho:  K-0  0.95F(1,8)  -  5.32 


Ho: K-0  ,  OK 


pnj  co  -4*  m 


Table  5.  Liquid  water  in  sn<x*pack  microwave-dilution  comparison, 

Analysis  of  Variance 
Field  Data 


Homogeneous  Snowpack 


Date 

X  (DL) 
(vol  Z) 

X  (FM) 
(vol  X) 

Deviation 

V 

v2 

Deviation 

v-v 

(v-v)2 

1 

1  May 

1.27 

3.30 

-2.03 

4.1209 

-1.67 

2.7889 

2 

7  May 

1.79 

1.70 

0.09 

0.0081 

0.45 

0.2025 

3 

8  May 

3.08 

4.20 

-1.12 

1.2544 

-0.76 

0.5776 

4 

10  May 

3.28 

4.15 

-0.87 

0.7569 

-0.51 

0.2601 

5 

10  May 

1.06 

1.02 

0.04 

0.0016 

0.4 

0.16 

6 

10  May 

2.27 

2.22 

0.05 

0.0025 

0.41 

0.1681 

7 

16  May 

3.19 

2.38 

0.79 

0.6241 

1.06 

1.194 

8 

16  May 

2.67 

3.4 

-0.73 

0.5329 

-0.37 

0.1369 

9 

23  May 

2.40 

3.02 

-0.6 

3.6 

-3.32 

0.1102 

0 

23  May 

4.23 

2.53 

1.7 

2.89 

1.97 

3.8730 

y 

-  -2.68 

0.00 

Model 

Residual 

DF 

Sum  of  Squares 
of  residuals 

lean 

square 

-  X  (FM) 

V 

-  X  (DL)  - 

V 

0  +  (e-6) 

10 

l  v2  -  10.551 

1.055 

-  X  (FM) 

V 

-  X  (DL)  =■ 

V 

K  +  (e-d) 

9 

l  (v-v)2  -  y.  13 

1.0936 

Difference 

1 

[A  -  0.718 

0.718 

F  »  0.657 


Ho:  K-0  0.95F(1,9)  -  5.12 


Ho:K-0  ,  OK 


csi  co  <r  m  \o 


Table  6.  Liquid  water  in  sncwpack  microwave-dilution  conqsarison, 

Analysis  of  Variance 


Laboratory  Data 


Homogeneous  Snowpack 


X  (DL) 
Date _ (vol  %) 


X  (FM) 
(vol  X) 


Deviation 

v  v2 


Deviation 
v-v _ 0 


24  Apr 

1.36 

10.0 

-8.64 

74.650 

0.112 

0.1254 

24  Apr 

1.24 

4.19 

-2.95 

8.7025 

5.802 

33.663 

25  Apr 

1.81 

5.55 

-3.74 

13.988 

5.012 

25.120 

1  May 

3.71 

29.55 

-25.79 

665.12 

-17.04 

290.29 

1  May 

4.64 

20.02 

-15.36 

235.93 

-6.608 

43.666 

1  May 

2.58 

13.02 

-10.42 

108.58 

-1.668 

2.7822 

7  May 

2.76 

8.75 

-5.99 

35.880 

2.762 

7.6286 

8  May 

3.62 

7.4 

-3.78 

14.288 

4.972 

24.721 

8  May 

2.68 

10.3 

-7.62 

58.064 

1.132 

1.2814 

10  May 

3.37 

6. 6 

-3.23 

10.433 

5.522 

30.492 

l 

*  -87.52 

0.00 

Residual 

Sum  of  Squares 

Mean 

Model 

DF 

of  residuals 

square 

(FM) 

-  Xv(DL)  - 

0  +  (e-6) 

10 

£  v2  ■  1225.6 

122.56 

(FM) 

-  Xv(DL)  - 

K  +  (e-6) 

9 

l  (v-v)2  -  459.66 

51.073 

Difference 

1 

-  765.98 

765.98 

F  -  15.00 


Ho:  K-0  0.95F(1 ,9)  -  5.12 


Ho:K-0  ,  Not  Valid 


Figure  Capitons 


Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 


1:  Comparison  of  electrical  path  length  model  with  data  of  Cumming 
(1952)  and  Sweeny  and  Colbeck  (1974)  (Colbeck,  1980). 

2:  Comparison  of  electrical  path  length  model  with  data  of  Ambach 
and  Denoth  (1972). 

3:  Comparison  of  electrical  path  length  model  with  Colbeck  (1980), 
Case  I  p3  »  0.32  Mg/tr  ew  ■  87.91. 

4:  Worst  case  percentage  error  of  microwave  measurement  of  liquid 
water  in  snowpack. 

5:  Insulated  snowbox  used  for  laboratory  intercomparisons. 

6:  Microwave  radiation  intensity  pattern  at  top  surface  of  the 
snowbox.  Relative  intensity  normalized  to  unity  at  center. 

7:  Example  of  liquid  water  distribution  in  the  snowbox;  measure¬ 
ments  are  in  percent  liquid  water  by  mass. 

8:  Microwave  response  to  snowpack  for  2-5  GHz  and  5-8  GHz  frequency 
bands. 

9:  Photograph  of  snowpit  used  for  field  measurements.  Note  several 
layers  on  the  pit  wall. 

10:  Pit  profile  for  April  19,  Runs  1,  2  and  3  field  measurements 
for  inhomogeneous  pack  are  depicted  on  the  right. 

11:  Pit  profile  for  May  1,  runs  4,  5  and  6  field  measurement  for 
inhomogenous  pack  are  depicted  on  the  right. 


Figure  12:  Photograph  of  snow  surface  after  being  sprayed  with  a  fine  mist 
of  dye.  Dark  splotches  are  areas  of  melting  snow. 
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